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Abstract 

We investigate the NLO QCD and the CP-violation effects in ttZ^ production at the 
Large Hadron Collider (LHC) in the minimal supersymmetric standard model(MSSM). 
Our calculation shows that the total NLO QCD correction in the framework of the CP- 
conserving MSSM significantly improves the scale uncertainty at the leading order, and 
the contribution from the pure supersymmetric QCD (pSQCD) correction can exceed 
-8% with the restrictions of 90 GeV < pip < 120 GeV and 120 GeV < p| < 150 GeV, 
where pip and are the transverse momenta of the top-quark and gauge boson, 
respectively. Our numerical results demonstrate that the pure supersymmetric QCD 
correction generally suppresses the total SM-like QCD correction in the CP-conserving 
MSSM, and tends to be a constant when either ti or g is heavy enough. We find also that 
the CP-odd asymmetry can reach 2.17 x 10~^, if the CP-phase angle really exists in 
the coupling of gluino-stop-top. 



PACS: 12.60.Jv, 14.70.Hp, 14.65.Ha, 12.38.Bx 
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I. Introduction 



Although the standard model(SM)[T],[2] has achieved great success in describing all the available 
experimental data, it suffers from some conceptional difficulties. That has triggered an intense 
activity in developing extension models. The supersymmetric (SUSY) extensions [3l IH |5l [6l [7]) 
rank among the most promising and well-explored scenarios for new physics at the TeV scale. 
Apart from predicting a light Higgs boson and describing the low energy experimental data very 
well, the SUSY models are able to solve various theoretical problems, e.g., the SUSY models 
may provide an elegant way to construct the huge hierarchy between the electroweak symmetry 
breaking(EWSB) scale and the grand unification scale. At present the minimal supersymmetric 
standard model(MSSM)[8] is regarded as the simplest and the most attractive one in the SUSY 
models. 

The direct evidence for the top-quark was presented in 1995 by the CDF and DO collabo- 
rations of the Fermilab Tevatron|9l |10]. This is considered to be a remarkable success of the 
SM. From that time on the top physics program has been turned to precise investigation for its 
properties. Since the top quark is the heaviest particle in the SM detected until now, it plays 
a special role in the mechanism of the EWSB, and the new physics connected to the EWSB 
may be found firstly through precise study of top-quark observables. The high accumulated 
top-quark events at the CERN LHC will open a new, rich field of top-quark phenomenol- 
ogy. Deviations of experimental measurements from the SM predictions, would indicate new 
non-standard top production or decay mechanisms. Therefore, the precise study of the top 
properties is one of the urgent priorities of the high energy experimental program. 

Beside the SUSY particle direct production, virtual effects of SUSY particles may induce 
deviations on observables from the SM predictions. If SUSY particles are really detected at 
the LHC, the comparison of precisely measured top-quark observables with the theoretical 
predictions including SUSY loop effects may yield additional information about the underlying 
model. Therefore, probing precisely the properties of the top-quark is an important goal at 
the LHC. In order to study precisely the top-quark physics within the SM and beyond at the 
LHC, it is necessary to give the theoretical predictions for top-quark observables including 
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higher-order corrections. In Ref . |1 1] S. Berge et ah, provide the predictions including the NLO 
SUSY QCD effects in the MSSM for the total production rate and kinematic distributions of 
polarized and unpolarized top-quark pair production at the Tevatron and the LHC 

Probing the couplings between the top quark and gauge bosons is another way to discover 
new physics. Until now there have been many works which devote to the observables related 
to the top-quark gauge couplings in the SM and beyond. The theoretical study of the effects 
of the top-quark and Z° gauge boson coupling at colliders was widely carried out. The cal- 
culations for the process e^e~ — t- tiZ^ at the leading-order(LO) and including next-to-leading 
order(NLO) QCD, electroweak corrections in the context of the SM were presented in Refs. 
[121 13], respectively, while CP- violating effects in e^e^ — ?■ tiZ^ process were studied in the 
framework of the two Higgs doublets model(THDM)[14j and with model independent effec- 
tive Lagrangian[T3]. The 77 — ?■ ttZ^ production channel has an outstanding advantage over 
e^e~ — )■ tiZ^ process in measuring tiZ^ coupling at the ILC, due to its relatively larger pro- 
duction rate[TSl [IT]. The NLO SUSY QCD corrections to the 77 — tiZ^ process at linear 
colliders, and the NLO SM QCD corrections to the tiZ^ production at the LHC are studied in 
Refs. [IS] and [IS], respectively. 

In this work, we calculate the production of the top-quark pair associated with a Z^ boson at 
the CERN LHC in both the leading-order (LO) and NLO QCD approximations in the framework 
of the MSSM with CP-conservation or CP-violation, and investigate the possible CP-violating 
effects contributed by the CP-phase in the couplings of gluino-stop-top predicted by the CP- 
violating MSSM. The paper is organized as follows: The description of the related theory about 
the CP-conserving and CP-violating MSSM is presented in section 2. In section 3 we outline the 
technical details of the related LO and NLO QCD calculations. In Sec. 4 we give some numerical 
results and discussions about the NLO SUSY QCD corrections and the possible CP-odd effect. 
Finally, a short summary is given. 
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II. Related theory of the CP- violating MSSM 

In the MSSM, each quark has two scalar partners called squarks, qi and QrIot qi and ^2)- The 
mass term of scalar quarks can be written as [20] 



->Cr^ = (^1 ^'r)MI[ ), (2.1) 



where A1| is the mass squared matrix of ql and qr, expressed as 

= ( "^S ""'1' 1 • (2.2) 
The diagonal and nondiagonal elements of this mass squared matrix are 



= M5 + m2 + m|(J^3-g,4)cos2/3, 



"^L = Mu,D + ^1 + Qqm^z^w COS 2/3, 



aqTUq 



rUg (Ag - ijru,D) , (2.3) 



where ruq, Qq and are the mass, electric charge and the third component of the weak isospin 
of the quark q, /i is the Higgsino mass parameter, Ag (g = u,d,c, s,t,b) are the supersymme- 
try soft-breaking trilinear coupling constants, Mg, M^ and M|, are the supersymmetry soft- 
breaking mass parameters of the left- and right-handed scalar quarks, and ru = I/td = cot f3 
are for the up- and down- type squarks, respectively. The mass squared matrix A^? can be 
diagonalized by introducing a unitary matrix TZ^. The mass eigenstates qi, q2 are defined as 



Then the mass term of scalar quarks can be expressed as 



<g 2 / qi 

where 



- = ( i ql ) ( I ) , (2.5) 



Mi' = WMpl^ t = ( < ^0^^ ) . (2.6) 



It is well known that the unitary matrix TZ^ can be parameterized as 



^q ^, cosOgC sin^ge*<^« , 

^ -sin^qC"^'^' cos^ge*<^« ^ ' 



where 6q is called as the mixing angle between the left- and right-handed squarks, and 20g is 
the phase angle of defined clS Clq — |Cig |e The masses of the squark mass eigenstates and 
the mixing angles acquire the forms as 

tan2^, = 4^^^^^. (2.8) 

Because of the large masses of the third generation quarks, the mixing effects of the third 
generation squarks are more significant than the first two generations. If we take the stop 
masses(mj^, m^^) ^^^d the stop mixing angle(6't) as the input parameters for the stop sector, the 
values of m^^, m^^ and \at\ can be obtained by adopting Eq. fl2.9p . 

2 2 /) 2 I • 2 /) 2 

m^^ = COS Bqm^^ + sm dqUi^^ , 

2 • 2 /I 2 I 2 /I 2 

m^-^ = sm 6'gmg^+cos fe'^mg-^, 
mg|ag| = sin cos (m?^ - m?J . (2.9) 

In the CP-violating MSSM, the SUSY soft-breaking trilinear coupling Aq and the Higgsino 
mass parameter /x can be complex. That makes having complex value. By using the param- 
eterization of the unitary matrix (Eq. (12.71) ). we obtain the squark current eigenstates (g^, 
in terms of the mass eigenstates (gi, q^) as 

gL = (gi cos^g - g2 sin^g)e^'^', q^ = (gi sin^^ + ga cos ^g)e-^'^'. (2.10) 

Normally the CP-violating effects in the MSSM from the gluino- squark- quark interactions 
are much more important than from the chargino and neutralino sectors due to the strong inter- 
action. We consider only the CP- violating effects induced by the ^ — ti.2 —t strong interactions. 
The Lagrangian for the gluino-stop-top couphngs is given by 

8 

a=l 
8 3 

= V2gsY,Y. ^-^-P {^imPL - t*hpPB)r + h.c, (2.11) 

a=l a,p=l 
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where Qs is the strong couphng constant, T° = {T^p) (a = 1, 8) are the 5'[/(3) generators, a, 
a, (5 are the color indices of gluino, top and stop separately, Pl,r = (1 =F7^)/2 and e = e"*'^^'^*^). 
Here 2(j)su(3) is the phase angle of Msu(3), the supersymmetry soft-breaking SU{3) gaugino 
(gluino) mass parameter, which is defined as 

Msuis) = |M5t/(3)|e'*'^^^(^). (2.12) 

By inserting Eg. (12 .101) into Eg. (12. lip , the mixing angle 6t and phase angle (pt may enter 
into the couplings, and the Lagrangian is expressed in terms of the stop mass eigenstates (^1,2) 
instead of the current eigenstates (tL,ij) as 

8 

Cg_i_t = \/2(7,^tT'^ [(tisin^i + t2Cos^t)e-<'^'+'^s^(3))p^ (2.13) 

a=l 

- (ti COS Ot - k sin Ot^ e^('^*+'^^^(^))PR g" + h.c. . 

As shown in this Lagrangian, only the combination of the phase angles 0t and (f)su{3) enters 
into the gluino-stop-top couplings. Therefore, we redefine this combination as 4>t, 

<Pt + (psuis) ^ (pt, (2.14) 
and obtain the conventional expression of the gluino-stop-top interaction Lagrangian as 

8 

Cg_i_, = V2g, J2 [ {ii sin Ot + ^2 cos Ot) e"*"^* Pl 

- (ti cos 9t - h sin 6^ c'^^'Pr g" + h.c. 
E E K^i^^ + V-ltPR)rU + h.c, (2.15) 



a=l 



2 8 

-I 

1=1 a=l 



where V,^^^ = iy/2g,T'' smete-''t'\ V-|^ = iv^^.T'^ cos ^fC"*'^*, V-|^ = -zv^c/.T'^ cos^te*"^* and 
VJ^^^ = iy/^gsT"" sin Ote'^'^* . There are similar expressions for other g — Qi — q{q = u, d, c, s,b) 
couplings involving CP-phase angles. Because in this work we consider only the CP-phase 
effects from the g — ii 2 — t couplings, we take 0t 7^ and 0^ = for g = u, d, c, s, b. 

In order to describe the CP-violating effects on the process, we take a definition of a CP- 
odd observable for the LHC, which is constructed to describe the distribution asymmetry of 
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(1) (2) (3) (4) 

Figure 1: The LO Feynman diagrams for the qq ttZ^ [q = u, d, s, c) partonic process. 

the azimuthal angle $ between pj and pj in the range of —180° < $ < 180°, i.e., 

$ = sgn [{pt - pt) ■ z] sgn [{pt X pi) ■ z] cos"^(pf ■ ), (2.16) 

where $ in Eq. fl2.16p comes from the modified definition of Eq.(14) in Ref.|21], and z is a unit 
vector of the z-axis direction along one of the incoming proton. The CP-asymmetry of angle $ 
is defined as 

^ Aa$ _ (t(180° > $ > 0°) - cr(0° > $ > -180°) 
* " ~W ~ cr(180° > $ > 0°) + (t(0° > $ > -180°) ' ^ ' 



The significance is defined as 



S=\^. (2.18) 



Then the CP-asymmetry effect may become observable at the Sa significance, if the integrated 
luminosity has a value larger than 

III. Calculations 

III..1 The LO cross sections for the partonic processes 

The contributions to the hadronic process of top-pair production associated with a boson 
at the LO, are from the partonic processes qq — )■ ttZ^ {q = u, d, c, s) and gg — t- ttZ^ channels. 
We use the 't Hooft-Feynman gauge in the following LO and NLO calculations. The LO 
Feynman diagrams for the subprocesses q{pi)q{p2) — ^ t{p3)i{pi)Z'^{p^), (g = u,d,c,s) and 
9{Pi)9{P2) t{P3)i{P'i)Z^{P5) in the MSSM are depicted in FigJTJand FigJ21 respectively. 
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(5) (6) (7) (8) 

Figure 2: The LO Feynman diagrams for the gg — > partonic process. 

The LO Feynman diagrams for all these subprocesses in the MSSM are the same as their 
corresponding ones in the SM. The exphcit expressions of the LO cross section for the partonic 
processes can be written in the form as: 



^ / dT,Y,\Mto? (3.1) 



^'^ 4|pi|\/^ 

where ij = qq, gg{q = u, d, c, s), the summation is taken over the spins and colors of initial and 
final states, pi is the c.m.s. momentum of one initial parton, and the bar over the summation 
recalls averaging over the spins and colors of initial partons. dr3 is the three-body phase space 
element expressed as 

dPa = {2nr5\p, +P. - R tMe- ^^''^ 

i=3 i=3 ^ ' * 

In the LO calculation for the parent process pp tiZ^+X we involve the contributions from 
partonic processes gg — > ttZ*^ and qq — > ttZ^ {q — u, d, c, s). Our numerical calculation shows 
the contribution to the LO integrated cross section from the partonic processes ss, cc — > ttZ^ 
is less than 3% at the LHC. Therefore, we consider only the NLO QCD corrections to the 
processes pp — > uu, dd, gg — >■ ttZ^ + X in the following NLO calculation. 

III.. 2 NLO QCD corrections to the partonic processes 

The NLO QCD correction in the MSSM(NLO SQCD) to each of the partonic subprocess qq 
tiZ^ {q = u, d) and gg ttZ^ consists of two independent parts. One is the so-called SM-like 



component, another is the pure SUSY QCD (pSQCD) component arising from the contributions 
of the virtual gluino one-loop diagrams. We adopt the dimensional regularization scheme in 
D = 4 — 2e dimensions to isolate the ultraviolet (UV) and infrared (IR) singularities. Then the 
total NLO SQCD corrections to the partonic subprocess qq — )■ tiZ^ {q = u, d) and gg — )■ tiZ^ 
, can be written as: 

^^Inlo = ^^SM-iike + ^^pSQCD^ i^J = ^d, gg). (3.3) 

In the MSSM, the so-called SM-like NLO QCD correction component is exactly equal to the 
NLO QCD correction in the SM, and we shall compare our results for the SM-like correction with 
those in Ref. [19]. The pSQCD correction component is UV and IR finite after renormalzation. 
The NLO SQCD correction includes the following contributions: 

► the virtual corrections to the partonic process qq{gg) — )■ ttZ^ . 

► the real gluon emission partonic process qq{gg) — )■ tiZ^g . 

► the real light- ( ant i) quark emission partonic process q{q)g — ?■ ttZ^q{q) . 

► the collinear counterterms of the PDF. 
(1) Virtual corrections in the MSSM 

In the MSSM, the virtual QCD 0{as) corrections come from the one-loop diagrams includ- 
ing self-energy, vertex, box and pentagon diagrams. In Figsj3]|6l we illustrate all the pentagon 
graphs for the partonic processes qq — )■ tiZ^ and gg — t- ttZ^ , separately. We take the defini- 
tions of the scalar and tensor two-, three-, four- and five-point integral functions presented in 
Ref.|23]. We use Passarino-Veltman method [2l] to reduce the N-point(A^ < 5) tensor functions 
to scalar integrals, and manipulate the 75 matrix in D-dimensions by employing a naive scheme 
as presented in Ref . [25] , which keeps an anticommuting 75 in all dimensions. The one-loop Feyn- 
man diagrams and the corresponding amplitudes are created by using FeynArts3.2 package [26]. 
and the scalar integrals are evaluate mainly by adopting the LoopTools-2.1 package [27] 128] . In 
order to cancel the UV divergences from both the SM-like and pSQCD one-loop diagrams, we 
should introduce some suitable counterterms. 
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Figure 3: The SM-like pentagon Feynman diagrams for the qq ttZ^ partonic process. 




Figure 4: The pSQCD pentagon Feynman diagrams for the qq — t- ttZ^ partonic process with 
the upper indices in t"'* and g*'* running from 1 to 2 respectively. 



9s gs + Sgs = gs + Sg 
tL ^ {l + ^6Zi)tL : 



SM-like 



+ 5g- 



pSQCD 



1 



1 

2' 

1 
1 

2 



1 I ^ f <: ryt,SM~Uke . <:ryt,pSQCD\ 



(3.4) 



where R and are the wave functions of top-quark and gluon, respectively. 

The counterterms of top-quark, gluon fields and top-quark mass are fixed by using on- 
mass-shell renormalization conditions [23j. For the renomalization of the QCD strong coupling 



constant Qs-, we use the MS scheme except that the divergences associated with top quark and 
the colored SUSY particle loops are subtracted at zero momentum[22]. The counterterm of the 
strong coupling constant includes the SM-like and pSQCD terms, which can be expressed as 
follows, 



5g. 



SM-like 



as(/ir 

An 



r nSM-like 



/3o 



1 1 ^ mi 
- -h - In ^ 
e 3 /i; 



(3.5) 
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(9b) ^ (10b) Z (lib) ^ (12b) Z 



Figure 6: The pSQCD pentagon Feynman diagrams for the gg — )■ partonic process. The 
upper indices in t^'* run from 1 to 2 respectively. The diagrams obtained by exchanging initial 
gluons are not shown. 
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An 



npSQCD ^ 2 i=l,2 ^ j=l,2 ^2 ^ 

+ ^ln^+ y -ln^+ y -In^ 

9? ;;2 19 /;2 19 /;2 



(3.6) 

where 

Pr-''-' = yAT, - In,, - I = -'-N^ - + 1), (3.7) 

with Nc = 3, nif = 5 hght flavors and 4 = + ln(47r) — Je- With the introduction of 
the CP-violating phase (pt in the MSSM, the renormahzed one-particle irreducible two-point 
functions for top-quark and gluon containing the contributions from pSQCD are deflned as 
follows [221 1 



It should be mentioned here that in the flrst equation of Eqs.f l3.8( ) the upper conjugation symbol 
I acts only on the CP-violating phase. The SM-like components for the top quark, gluon self- 
energies, wave function and top mass counterterms, can be found in many references, such as 
Ref . |23] . Here we present only the related pSQCD component expressions for unrenormalized 
top quark, gluon self-energies and counterterms. 

^fsQCDip") = -^,9l (cos^ Oi B,[p\ ml m|] + sin^ 9-, B,[p\ mj, m|]) , (3.9) 

^pIqcd{p') = -%9l (sin' 9-, B,\p\ m|, m|] + cos^ 9-, B,\p\ m|, m|]) , (3.10) 

^fslcDip") = -%9lmg (sin% cos^,- e'^^*) (5o[p^m?,m?J - So[p^m^,m|]) , (3.11) 

^%Up') = -^29>9 (sin^f cos^, e^^*) (5o[/,m|,m|] - i?o[p^m|,m|]) , (3.12) 
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Since we don't need the longitudinal part of gluon self-energy in the following calculation, 
we give the explicit expression for ^^pSQcnip"^) only. 



^pSQCdKP ) 



9t 



167r2 

2 s,t,b 

+ E E (^K]-25oo[p^ 

i=l q=u,d,c 



{3D [(2 - D)Boo + mlBo - p^Bn + B,)] [p^, ml 

])}, (3.13) 



2 2 2 



where D = 4 — 2e and the definitions of the two-point integrals are adopted from Ref . |23] . In 
the SU{3) group, Cp = {N^ — l)/{2Nc). By using the relevant on-mass-shell renormalization 
conditions and imposing the real condition on the right-handed top-quark field renormalization 
constant, 6Zj^ = 5Z^[29], we obtain 



6m 



pSQCD 

t 



^tRe^%CDi'rn^t) + mtRe^'-p^QCDimt) + ReJ:''pSQCDii^t) + R^^pSQCD 



(t)S 



m 



t)i 



^ ^t,pSQCD 
L 



^pSQCoij^^t) 



mt 



<t)L 
pS( 

_d_ 

.{t)S 



R'^'^plqcoi'^l) ~ R^'^isQCoi^t , 

(t)R 



+ ReT.';itQCD{p') + Re^fsUoiP 
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(3.14) 



(3.15) 



^ 2^t,pSQCD 
R 



~ ^pSQCDV'h . 



{t)S t 
pSQCD 



^^t^fsQCDip') + ^t^sZcD^P 



(t)R 



(3.16) 



^ ^pSQCD 



~ d^%cnip') 



-Re 



(3.17) 



where Re only takes the real part of the loop integral functions appearing in the self-energies. 
The renormalized amplitudes of all the NLO QCD virtual corrections to the partonic processes 
qq — )■ tiZ^ and gg — t- tiZ^ in the MSSM are expressed as 



AMIL = AMl,f + AMIU + ^Mll + AM^^, + A>1^L„„ (tj = uu, dd, gg), (3.11 



box 



'■pent 



14 



where AMl^^ip AMU^, AMl^^^, AMp^nt ^-^^count represent the amphtudes for self-energy, 
triangle, box, pentagon and counterterm diagrams, respectively. Then we can get the UV-finite 
virtual NLO QCD correction component Aa*-'-^ as 

Aal = — ^ / dr3^Re(Al^o x A>ll). (3.19) 

The definitions of the notations appeared in above equation are the same with those in Eq. fl3.ip . 
(2) Real gluon and light- (anti) quark emission corrections 

In the MSSM, the real 0{a^ correction processes involve the real gluon emission and real 
light- (anti) quark emission processes which are listed as follows: 

g(Pi) + g(P2) ^ t{vz) + t(p4) + Z\v^) + g{v^) (3.20) 
9{Vx) + 9{V2) ^ t{vz) + t(P4) + Z\v^) + 9{V^) (3.21) 

g(g)(Pi) + 9{V2) -> t(P3) + t(P4) + Z\v^) + g(g)(P6). (3.22) 

Because of the IR singularities involved in these processes, we use the two cutoff phase space 
slicing method (TCPSS) to perform the integration over the phase space of these real emission 
processes. [31] . In our calculations, the real gluon emission correction to each of the processes 
ij — )■ [ij = uu, dd, gg) contains both soft and collinear IR singularities, which are involved 

in soft gluon region(i?6 < 5sV^/2) and hard gluon region(i?e > 6s\^/2) respectively. The hard 
gluon region is also divided into the hard collinear region (HC) and the hard noncoUinear 
region (HC) with < 5c and > 6c {pi are the momenta for q and q). Each of the 

real light- (anti) quark emission processes contains only collinear IR singularity, and can be dealt 
with in the hard collinear region(HC) too. In the HC region, the real emission corrections, 
Aa^, where kl = qq, gg,qg,qg, {q = u,d), are finite and can be calculated numerically with 
general Monte Carlo method. After summing the virtual and real gluon/ (anti) quark radiation 
corrections, the remained collinear divergence can be cancelled by that in the NLO PDFs. Then 
the finite total NLO QCD correction to the pp — )■ + X process can be obtained. 
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III. .3 Total cross sections for the pp — > ttZ^ + X process 

The LO, NLO SQCD corrected hadronic cross sections for pp tiZ^ + X in the MSSM can 
be written as: 

cc,ss,gg ^ i 

daLo,SNLo{pp -^ttZ° + X) = / dxidx2 

x[Gi(xi,/i/)Gj(x2,/i/)(ia^^o,57VLo(^i^2v^'/"r) + (1 ^(2)p3) 

We adopt the CTEQ6L1 and CTEQ6m PDFs[32j in the LO and NLO calculations respectively, 
except in the specific case for numerical comparison. The (3"^o SA^Lol^i^aV^, yUr) are the LO, 
NLO SQCD corrected cross sections with the partonic colliding energy = X\X2\fs for the 
partonic processes of ij — ?■ ttZ^{ij = gg,qq) in the MSSM. Throughout our evaluation, we 
equate the factorization and renormalization scales and define fi = fij = fir- 

In the MSSM, the NLO SQCD corrected partonic cross sections can be expressed as below: 

dcrgj^^Q{xiX2y/s, /i) = (ia^Q(xiX2A/s, /i) + dAa''^,^^Q{xiX2y/s, fi) (3.24) 

where ij = gg,qq, A(3"^^^q(xi, X2, /i) denotes the total NLO QCD correction in the MSSM 
to the corresponding LO partonic cross section. In this work we ignore reasonably the NLO 
QCD corrections to the partonic processes cc, ss — )■ tiZ^ due to the luminosity suppression, 
i.e., A(T^'^^^Q(xiX2A/i, /i) = 0. Then the full NLO QCD corrections in the MSSM to the process 
pp — > tiZ^ + X at the LHC can be divided as: 

Ao'SNLO = A(^SM-like + A<^pSQCD- (3.25) 

The later part arises from the virtual correction of the diagrams with gluino/squark in loops. 
Note that the cross sections for the parent process pp — > ttZ^ + X asu-iike should be the same 
as the corresponding ones, asm, in the SM[T9l [33] . 

IV. Numerical Results and discussions 

In this section, we present numerical results of the NLO QCD corrections to the process pp — )• 
ttZ^+X in the MSSM at the LHC. The numerical results for the LO and the NLO SM-like QCD 
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/i 


our (TLo(pb) 


our (TAfLo(pb) 


o"LO (pb) in Ref.|33J 


o'NLoip^) ill Ref. [33] 


/io/4 


1.0779(8) 


1.216(5) 


1.078 


1.213 


/io/2 


0.8083(6) 


1.095(4) 


0.808 


1.093 




0.6198(3) 


0.975(4) 


0.620 


0.973 



Table 1: The comparison of our numerical results for the pp — )■ tiZ'^ + X process in 
the SM at the LHC with those in Ref. [33], where the LO and NLO QCD corrected 
cross sections for different energy scale values(/i = fij = /i^) are listed with the 
relevant parameters and the PDFs being the same as in Ref. [33], i.e., fio = 2mt+mz, 
mt = 170.9 GeV, mz = 91.19 GeV, mw = 80.45 GeV and the MRST PDFs. Note 
that the numerical results are contributed by the subprocesses qq, gg — ?■ tiZ^ with 
q = u,d,c, s at the LO, but q = u,d at the NLO. 

corrections have been compared with the data presented in Tabel 1 of Ref. [33] . Both result sets 
are listed in Table [H There we employ the MRST PDFs[34J and the input parameters which 
were used in Ref. [33]. The agreement between them can be seen obviously from the table. 

In the following numerical calculations in the frameworks of the SM and MSSM, we define 
fiQ = 1711 + iTLz/'2 and take CTEQ6L1 PDFs with an one- loop running as in the LO calculation 
and CTEQ6M PDFs with a two-loop as in the NLO calculation [35]. The number of active 
flavors is Nf = 5 and the QCD parameters are A^^ = 166 MeV and = 227 MeV for the 
LO and NLO calculations, respectively[35]. We ignore the masses of u-, d-, and s-quarks in our 
calculations. The other SM parameters are taken as|35j. 

aewimz^) = 1/127.918, = 171.2 GeV, nib = 4.2 GeV, = 1.3 GeV, 
mz = 91.1876 GeV, mw = 80.398 GeV, sin^ Ow = I - m^wlm\ = 0.222646. (4.1) 

For the SUSY parameters of the scalar top sector in the CP-conserving MSSM, we take 
the top squark masses(m(-^, m^J and their mixing angle(6'i) as input parameters, and adopt 
Eq. 02.91) to calculate the m?^, fn'j^, at, and sequentially use Eq. fl2.3p to get the SUSY soft- 
breaking parameters Mq^ and Mj/^ for the stop sector by assuming tan/3 = 10 (Here we take 
tan /3 = 10 arbitrarily, because the NLO SQCD cross sections for pp — )■ ttZ^ + X do not directly 
related to tan/3, and in the range of 1 < tan/3 < 50 the dependences of m^^ and m^^ on tan/3 
are rather weakfTTj). Since the gauge invariance in the MSSM requires Mq^ = Mq^, we need 
only to fix two additional parameters for the scalar bottom sector. We neglect the mixing for 
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the scalar bottom sector and assume m^^ = fni^, the masses of 61^2 can be obtained by the 
following equation, 



m 



+mi + ml{ll - Q,si,) cos 2/3. (4.2) 



bL 



Furthermore, we take rrig = 200 Gel/[3S], 4>g = (for the u-, d-, c-, s-, b-quarks) and let 
CP-phase = 0^ 7^ as a free parameter in the CP-violating MSSM. For the first two squark 
generations, we take the assumption of the universal squark mass as used in Ref. [TT] i.e., 

rriu^ = rriu^ = m^^ = mj^ = mg^ = = mg^ = m^^ = nig, (4.3) 

and set rriq = 1 TeV in the following calculations. 

To make the demonstration of the correctness of our calculation for the integrations over the 
phase space of the 4-body final-state real emission processes, we have checked the independence 
of the SM-like NLO QCD correction component Ao^^J^q^" of the process pp uu tiZ^ + 
X on the two cutoffs 6s and 6c separately shown in FigsJ71^a-b) and FigsJHl^a-b) separately. In 
identifying the Aa^^^p'^'^ independence on (5^, we fix (5c = 1 x 10~^ and vary 6s from 4 x 10"'* 
to 8 X 10~^. For probing the Aag^J^Q^'' independence on 6c, we take 5^ = 1 x 10~^ and let 
6c running from 1 x 10^^ to 4 x 10~^. These four figures show that although the three-body 
correction(Acr^^)) and four-body correction(Acr*^^)) are strongly related to the cutoffs 6s and 6c-, 
the final total SM-like NLO QCD correction Aa^^^^'"' to the pp uu tiZ^ + X process, 
which is the summation of the three-body term and four-body term, i.e., Aa^^^Q^^ = Acr*-^^ + 
Aa^^\ is indeed independent of the cutoffs {6s and 6c) within the statistical errors. 

The LO and NLO QCD corrected cross sections in the CP-conserving MSSM for the process 
pp — )■ tiZ^+X as the functions of the renormalization and fatorization scales in the 
way at the LHC, are demonstrated in FiglHl^a). The curve for the component of the SM-like 
NLO QCD corrected cross section is also shown there for the comparison. The corresponding 
K-factors for the NLO SQCD and SM-like QCD corrections, which are defined as K{fi) = 
(^SNLoil^)/crLoilJ^) and KsM-nkeifJ') = (^sN£o'"'il^)/^Lo{l^), are drawn in FigJHKb), respectively. 
There we take fiQ = nit + rrig = 200GeV and /i running from /io/5 to 3/io- As we 

know, the scale dependence of the PDFs for the incoming u- and rf-quarks are significant, 
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Figure 7: (a) The dependence of QCD NLO correction parts to the pp ^ uu ^ tiZ^+X process 
on the soft cutoff ds at the LHC with the colhnear cutoff 6c = Ix 10^^ and ^/s = 14 TeV. (b) The 
amplified curve for the total QCD correction ^cTsnlo^'^ the process pp ^ uu ^ tiZ^ + X , 
where it includes the calculation errors. 
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Figure 8: (a) The dependence of the QCD NLO correction parts to the pp ^ uu ^ ttZ° + 
X process on the colhnear cutoff Sc at the LHC with Ss = 1 x and ^/s — 14 TeV. (b) The 
amplified curve for the total QCD correction ^cfsnlo^'^ process pp ^ uu ^ ttZ^ + X , 

where it includes the calculation errors. 
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even up to the NLO the scale dependence of the results are not too small. Since the ttZ'^ 
production at the LHC is a QCD process at the LO, the scale dependence of the NLO SM- 
like QCD corrected cross section (dotted curve in Figj9]^a)) is less than that of the LO cross 
section for the pp — ttZ^ + X process. On the other hand, the pSQCD correction is to some 
extend a LO contribution, because this process does not involve the supersymmetric strong 
coupling Qs at the LO. Therefore, the pSQCD correction induces some more scale dependence 
to the pp — > ttZ^ + X process. Considering the fact that pSQCD correction is quite small 
comparing with the SM-like QCD correction demonstrated in FigJHI^a), the scale dependence 
of the full SQCD corrected cross section is similar to that of the SM-like QCD corrected cross 
section, which is much less than that of the LO cross section. Actually, FiglHI^a) demonstrates 
that the LO cross section is strongly correlated with the energy scale fi, while the NLO QCD 
corrections obviously improve the scale uncertainties in both the CP-conserving MSSM and the 
SM. Comparing with the SM-like NLO QCD correction, it can be seen the pSQCD corrections 
cancellate the correction part from the SM-like QCD, and the NLO pSQCD correction to the 
total cross section can exceed —4.75% in our chosen parameters space. FiglHIa) shows that the 
total NLO SQCD K-factor changes from 0.48 to 1.63 as the scale running from /xo/5 to S/iq. 
In the following calculations we set fi = fiQ. 

Fig JTUT a) shows the LO and NLO QCD corrected cross sections in the CP-conserving MSSM 
and the SM for the process pp — ?■ ttZ^ + X as the functions of the gluino mass {rrig) at the 
LHC by taking m^-^ = lOOGeV", m^-^ = GOOGeV, and the mixing angle dt = —tt/A for the stop 
sector. The corresponding K-factors for the NLO SQCD and SM-like QCD corrections versus 
rrig are depicted in FigJTUlfb). respectively. Figs ITUT a) and [TUT b) demonstrate that although 
the SM-like curves have no relation with gluon mass, the curves for NLO SQCD corrections 
is clearly related with rrig in the region of rrig < 300 GeV. The NLO SQCD corrected cross 
section (and K-factors) approaches a constant when rUg > 400 GeV due to the decoupling effect 
for heavy gluino exchanging. We can see that the NLO pSQCD correction is non-zero when 
rrig > 400 GeV, because of the relatively light mass of ti in loops (m^-^ = 100 GeV). FigJTUlfb) 
shows when we take rUg = 100 GeV, the pSQCD relative correction can reach to —8.56% and 
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Figure 9: (a) The dependence of the LO and NLO cross sections in the CP-conserving MSSM 
and the SM on the energy scale at the LHC. (b) The total NLO QCD K-factors for the 
process{K sNLoifi') = o'SNLo{fJ')/<^Lo{fi')) and the NLO QCD K-factors of the SM-like part 
{KsM-iikeifJ') = o'sNLo'^^if^) /'^LoifJ')) for the pp —J- + X process versus the energy scale at 
the LHC. 

K-factor of the total SQCD correction will be 1.281, while when we fix rrig = 200 GeV, we 
get —4.75% for the corresponding pSQCD relative correction and 1.319 for the K-factor of the 
SQCD correction. 

We present Fig JTlT a) to show the relations between the LO and NLO QCD corrected cross 
sections in the CP-conserving MSSM and SM for the process pp — )■ tiZ^ + X at the LHC as 
the functions of mj^, one of the input parameters for stop sector. Here we take rrig = 200 GeV, 
the other two input parameters for stop sector are set as m^^ = QOOGeV, 6t = — vr/4, and 
the other SUSY parameters are obtained as explained above. The corresponding K-factors 
for the NLO SQCD and SM-like NLO QCD corrections versus are depicted in Fig JTTT b). 
respectively. Again we see from Figs JTlT a) and (b) that the SM-like curves do not show the 
dependence on m^^, but the NLO QCD corrections in the MSSM obviously rely on m^^ in the 
region of m^^ < 400 GeV. While in the region of m^^ > 400 GeV the NLO QCD corrections and 
corresponding K-factors in the MSSM tend to be constant respectively, because there exists 
the decoupling effect of heavy ti in loops, and the nonzero NLO pSQCD correction is induced 
by the relatively light mass of gluon in loops(mg = 200 GeV). From Fig JTTT b) we can see that 
with the mass of ti running from 100 GeV to 500 GeV, the NLO pSQCD relative correction 
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Figure 10: (a) The LO and NLO QCD corrected cross sections in the CP-conserving MSSM 
as the functions of the gluino mass rrig at the LHC (b) The corresponding total NLO QCD 
K-factor in the MSSM for the process pp tiZ^ + X {KsNLo{mg) = ^^j^l^) and the 
NLO QCD K-factor for the SM-hke part for the pp — ?■ ttZ^ +X process versus the gluino mass 
at the LHC. 

varies from —4.75% to —0.25%, and the K-factor of the total NLO SQCD changes from 1.319 
to 1.364. 

The LO and NLO QCD corrected differential cross sections of the transverse momenta 
for top quark and boson in the CP-conserving MSSM and the SM for the process pp — )• 
ttZ^ + X at the LHC, are drawn in Fig JTW a) and Fig{T3]^a), separately. In these plots we take 
mi^ = lOOGeV^, m^^ = GOOGeV and the mixing angle 9t = — 7r/4 for stop sector, nig = 200 GeV, 
and the other SUSY parameters are set to have the values explained above. Their corresponding 
K-factors(K(pT) = "'l^y ° / ^re depicted in FigsHSl^b) and FigHSj^b), respectively. There 
we take m~g = 200 GeV. Both figures FiglEl^a) and FigHSj^a) show that the SM-like NLO 
QCD corrections enhance the differential cross sections of the transverse momenta for the top 
quark and Z° boson in whole plotted range, while the NLO pSQCD correction part suppresses 
the SM-like QCD correction slightly. We can obtain from Fig JT2l^ b) and Fig JTST b) that the 
relative corrections from the pSQCD can be —8.56% and —8.12% when pip ~ 105 GeV and 
Pp ~ 135 GeV respectively. And in these two figures there exist obvious distortions for the 
two NLO SQCD K-factor curves compared with the corresponding SM-like ones. Those curve 
distortions are caused by the resonant effect of gluon self-energy. 
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Figure 11: (a) The LO and NLO QCD corrected cross sections in the CP-conserving MSSM as 
the functions of the hght scalar top-quark mass m^^ at the LHC. (b) The corresponding total 
NLO QCD K-factor in the MSSM for the process pp ttZ^ + X (KimiJ = ^^f^/f^) and 

the NLO QCD K-factor for the SM-like part for the pp ttZ^ + X process versus m^^ at the 
LHC. 




Figure 12: The LO and NLO QCD corrected distributions of the transverse momentum of the 
top quark for the process pp — )> ttZ° + X in the SM and CP-conserving MSSM at the LHC and 
the corresponding K-factors(i^(p^) = ^^^r^/^^^) versus pip. (a) the differential cross section 
of the transverse momentum, (b) the corresponding K-factors. 
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Figure 13: The LO and NLO QCD corrected distributions of the transverse momentum of the 
Z°-boson for the process pp — )■ + X in the SM and CP-conserving MSSM at the LHC and 
the corresponding K-factors(-ft'(pf ) = ^^^^/^^) versus . (a) the differential cross section 
of the transverse momentum, (b) the corresponding K-factors. 

In FigJTH we define K{M^tt)) = ^^^m^Iim^ ^'^^ P^°^ curves for the differential cross 

(tt) (it) 

sections and corresponding K- factor as the functions of the top-pair invariant mass M(ti^^ where 
we take m^^ = lOOGel^, mi^ = GOOGeV, 9t = —tt/A for stop sector, and rrig = 200 GeV. There 
we see the K{M(^tt^) distribution demonstrates the characteristic effects, which are shown already 
on the K-factor curves in FigJTWb) and Fig JTST b). The distortion of the K-factor distribution 
curve for M(^tt) is located at the vicinity of M(^tt) ~ 2mg = 400 GeV where the K-factor of the 
NLO SQCD correction reaches the value of 1.457. It exhibits exactly that the resonance effect 
of the gluino-pair threshold in the gluon self-energy induces the curve distortion. 

If the CP-violating MSSM is true, the CP-odd effects for the process pp — t- ttZ^ + X at 
the LHC would be demonstrated through the cross section deviation from the CP-conserving 
MSSM, and shown from a nonzero CP-odd observable Aip defined in Eg. (12. 171) . We plot the 
NLO SQCD corrections to the cross sections of pp — tiZ^ + X process as the functions of 
the CP-phase (p FigJTSlfa). and the corresponding K-factors are drawn in Fig JTST b). In 
these two figures we set (p = (pt and (pq = 0, {q = u,d, s,c,b), rrig = 200 GeV, for the stop 
sector we take two set of input parameters: (1) {m^-^, m^-^, = {lOOGeV^, GOOGeV, — ^}(dash- 
dotted-curve); (2) {mt^,mi^,et} = {250^61/, SOOGeV', -f}(dotted-curve), and the other SUSY 
parameter values are explained above. The curve for the SM-like NLO QCD correction (fuU- 
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Figure 14: The LO and NLO QCD corrected differential cross sections of the top-pair invariant 
mass M(tf) and the corresponding K-factors in the CP-conserving MSSM at the LHC. 

hne) in Fig JTST a) shows the SM-hke correction part does not depend on the CP-phase, while 
the curves for NLO SQCD correction part vary as cosine wave of (j). The curves for K-factors 
of the SM-like QCD and the SQCD correction in Fig JTST b). show the similar behaviors with 
the corresponding ones in FiglTSlfa). 

We adopt the CP-asymmetry parameter (^$) definition in Eq. fl2.17p and depict Aq, as the 
function of CP-phase in FigJT6l where nig = 200 GeV and the solid curve is for taking 
{mi^,mi^,9t} = {lOOGeV, GOOGeV, -f }, the dash-dotted curve {m^^, m^-^, 0^} = {250GeV, SOOGeV, -f }. 
It shows that the two curves for Aip are sinusoidal with respect to $, and the absolute value 
of for {mj^, m^j, 9t} = {lOOGeV, GOOGeV, —j} can reach its maximum of 2.17 x 10"'^ when 
0= = 1,3,5,7), while it can be 1.57 x lO^^ for {mi^,mi^,9t} = {250GeV, SOOGel/, -f }. 

From Eg. (12. 191) we can see that if we assume the total cross section of the process pp — )■ 
tiZ^ + X is 1 pb and the integral luminosities are larger than 212 fb~^, 1911 fb~^ or 5309 fb~^ 
we may observe the CP-violating effect induced by at la, 3a or 5a significance for the case 
of {mi^,mi^,9t} = {lOOGeV, GOOGeV, -f }, and 406 fb-\ 3651 fb'^ or 10142 fb'^ for the case 
of {mt^,mi^,9t} = {250GeV, SOOGeV, -f }, respectively. 
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Figure 15: (a) The LO and NLO QCD corrected cross sections as the functions of the CP-phase 
<p{= 4>t) for the process pp tiZ^ + X at the LHC. (b) The corresponding total NLO QCD 
K-factors versus CP phase (f). 




Figure 16: The CP-violating asymmetry parameter A<s> as the functions of the CP-phase angle 
for the process pp — )■ tiZ^ + X at the LHC, where the solid curve is for {mf^,mf^,9t} — 
{IQOGeV, GOOGeV, -7r/4}, and the dashed curve for {mj^, m^^, 9t} = {250^6^, SOOGeV, -7r/4}. 
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V. Conclusion 



In this paper, we study the total NLO QCD corrections for the process pp — > ttZ^ -\- X m the 
MSSM at the LHC. Our results show that both the total NLO QCD corrections in the CP- 
conserving MSSM and the SM improve the LO scale uncertainty. We provide the NLO QCD 
corrected distributions of transverse momenta of the top-quark and boson at the LHC in 
the CP-conserving MSSM and the SM. There we can see that the total NLO QCD corrections 
can modify significantly the LO cross sections respectively. The pSQCD corrections to the 
process pp — >■ ttZ^ -\- X can be beyond —8% when we take rUg = 200 GeV and restrict the 
top-quark QOGeF < < 120 GeV or the Z° boson 120 GeV < p| < 150 GeV. And we see 
that the K-factor is sensitive to the value of rrig or mi^ in the relatively lighter mass region of g 
or ti- Furthermore, we find that if the CP- violating phase really exists in the scalar top mixing 
matrix or in the Majorana mass term of the gluino predicted by the CP-violating MSSM, the 
CP- violating effect, described by CP- violating asymmetry parameter A^, can be expected to 
be the order of 10~^ and reach the maximal value 2.17 x 10~^. Therefore, testing CP- violation 
induced by CP-phase in pp ^ tiZ^ + X process could be an interesting task at the LHC. 
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